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Abstract 



The reaction 7p — ► 7r°7r°p has been measured using the TAPS BaF2 calorimeter 
at the tagged photon facility of the Mainz Microtron accelerator. Chiral pertur- 
bation theory (ChPT) predicts that close to threshold this channel is significantly 
enhanced compared to double pion final states with charged pions. In contrast to 
other reaction channels, the lower order tree terms are strongly suppressed in 2n° 
photoproduction. The consequence is the dominance of pion loops in the 2tt° chan- 
nel close to threshold - a result that opens new prospects for the test of ChPT and 
in particular its inherent loop terms. The present measurement is the first which is 
sensitive enough for a conclusive comparison with the ChPT calculation and is in 
agreement with its prediction. The data also show good agreement with a calculation 
in the unitary chiral approach. 
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1 Introduction 



In the energy regime where excitations of the nucleon or properties of the 
lowest lying mesons are studied, the perturbative ansatz to solve QCD fails, 
because the strong coupling constant a s is too large. A different approach 
exploiting the approximate Goldstone boson nature of the pion has be en de- 
veloped, namely chiral perturbation theory (ChPT) ()Wei79l: lGas84| ). This 
effective field theory has been extended to the nucleon sector (HBChPT 1 ) 
f),Ten9lt lBer92h . Historically, it was a great success that ChPT could clarify 



the disagreement between the old low-energy theorems (LET s) in describin g 
the s-wave threshold amplitude E 0+ in the reaction jp —>■ 7i°p 



Additional contributions due to pion loops had to be added to the old LETs, 
resulting in a lowe r value for the E 0+ amplitude and agreement with mea- 
sured data ( Ber90j ). In gener al, ChP T is in good agreement with experiments 



describing n — N scattering (jFetOOh 



From the study of irir production processes, complementary information to 
the study of the single pion photoproduction channels can be gained. The 
extention of ChPT to nn photo- and electroproduction has led to the finding 
that the cross section for final stat es with two neutral pions is dramatically 



enhanced due to chiral (pion) loops (JBer94j) which appear in leading (non van- 
ishing) order g 3 . This is a counter-intuitive result, since in the case of single 
pion production the cross sections for charged pions are considerably larger 
than those with neutral pions in the final state. This situation is not changed 
when the ChPT calculation is extended by evaluating all nex t-to-leading or- 
der terms up to order M% in the threshold amplitude ( Ber96l ). Exploring this 



situatio n in more detail, the two pion channel exhibts the following properties 
f)Ber9fih : Born type contributions start at order M w and are very small. Tree 
diagrams up to order q 3 are zero due to threshold selection rules or pairwise 
cancellation. Only at order q 4 do tree terms proportional to the low energy 
constants Cj give a moderate contribution. In a microscopic picture these tree 
terms (a a) subsume all s-, t- and u-channel resonances in the 7r°p scattering 
amplitude (e.g. the A(1232)). The largest resonance contribution at order M% 
comes from the Pn(1440) resonance via the N*Nirir s-wave vertex. Possible 
double A graphs as well as loop diagrams containing a photon coupling to 
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a K + — S/A pair were estimated and found to be negligible. Fourth order 
loop diagrams (g 4 ) provide only a moderate contribution. All the coefficients 
of the resulting threshold amplitude were taken from the literature, tt — N 
scattering and, in case of the s-w aye Pi i (1440) to tttt coupling, from an anal- 
ysis of the reaction ttN — > ttttN ( Ber95| ). Adding all contributions together, 
the astonishing result is that the yield of the leading order loop diagrams (g 3 ) 
is approximately | of the total 2tt° strength. This fact makes this channel 
unique, because unlike in other channels where the loops are adding some 
contribution to the dominant tree graphs, here they dominate. Consequently 
the 2tt° channel provi des a y ery sensitive method to study these loop contri- 
butions to ChPT. In (|Ber96h . the following prediction for the near threshold 
cross section was given: 



/ e — E thr \ 2 

a tot (E„) = 0.6 nb ^ ?-) (1) 

K l! V 10 MeV J V ; 

where E y denotes the photon beam energy and E* is the production thresh- 
old of 308.8 MeV. Actually, the uncertainty of the coupling of the Pn(1440) 
to the s-wa ye tttt c hannel was a limiting factor for the accuracy of the ChPT 
calculation ( Ber96l ). For the most extreme case of this coupling, an upper limit 



for this cross section was deduced by increasing the constant in Eq. 1 from 0.6 
nb to 0.9 nb. 

To complete the overview of theoretical calculations of the reaction 7^ — > 
7r°7r°p close to threshold, it is noted that this chan nel is a lso described in 
a recent version of the Gomez Tejedor-Oset model (|Tei96h . This model is 



based on a set of tree level diagrams including pions, nucleons and nucleonic 
resonances. In a recent work, partic ular em phasis was put on the rescattering 
of pions in the isospin 1=0 channel (lRoc02h . Double pion photoproduction via 
the A Kroll-Rudermann term is not possible for the 2tt° final state. In the case 
of a tt~tt + Kroll-Rudermann term, the charged pions can rescatter into two 
neutral pions generating dynamically a tttt loop. This effect nearly doubles the 
cross section in the threshold region and is regarded by the authors as being 
reminiscent of the explicit chiral loop effect described above. Nevertheless, the 
cross section calculated with this model is significantly smaller than the ChPT 
prediction. 

In the past, two measurements of t he react i on 7p — > tt°tt°p below 450 MeV 
beam energy have been carried out (lHae97l : IWolOOi ). The second experiment 



showed an improvement in statistics by almost a factor 30. Nevertheless, in 
the threshold region the cross section still suffered from large statistical un- 
certainties (see Fig. 2). 

In contrast to previous analyses, we did not extract the cross section from 
events in which only three of the four decay photons were detected. Such 
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analyses are not kinematically overdetermined and, close to threshold in par- 
ticular, the extracted cross section can be slightly con tamina ted by the re- 
action 7p — ■> 7r°7p, which has recently been measured ( Kot02j ). This present 
measurement of the 2ir° photoproduction at threshold is the first for which 
comparison to theoretical calculations is conclusive. 



2 Experimental Setup and Data Analysis 



The reaction 7£> — > 7r°7r°p was me asured at the Mainz Microtron (MAMI) 
electron accelerator (|Wal90l : lAhrii ) using the Glasg ow tagged phot on facility 
( Ant9ll : lrh3J9fih and the photon spectrometer TAPS (lNov9lt lo"ab94[) . A quasi- 
monochromatic photon beam was produced by bremsstrahlung tagging, in 
which the photon energy is determined by the difference between the energies 
of the incident electron and the residual electron following bremsstrahlung be- 
cause the energy transfer to the atomic nuclei of the radiator foil is negligible. 
The photon energy covered the range 285-820 MeV with an average energy 
resolution of 2 MeV. The photon flux was of the order of 0.5 MHzMeV" 1 at 
photon energies of 300 MeV. The TAPS detector consisted of six blocks each 
with 62 hexagonally shaped BaF2 crystals arranged in an 8x8 matrix and 
a forward wall with 138 BaF2 crystals arranged in a 11x14 rectangle. Each 
crystal is 250 mm long with an inner diameter of 59 mm. The six blocks were 
located in a horizontal plane around the target at angles of ±54°, ±103° and 
±153° with respect to the beam axis. Their distance to the target was 55 cm 
and the distance of the forward wall was 60 cm. This setup covered ~40% of 
the full solid angle. All BaF2 modules were equipped with 5 mm thick scintil- 
lation plastic dE/dx detectors to allow the identification of charged particles. 
The liquid hydrogen target was 10 cm long with a dia meter o f 3 cm. Further 
details of the experimental setup can be found in ref. (|Kot03l ). 



The 7p — > ir°7r°p reaction channel was identified by measuring the 4-momenta 
of the two 7T° mesons, whereas the proton was not detected. For a three-body 
final state this provides kinematical overdetermination and hence an unam- 
biguous identification of this reaction channel. The tt° mesons were detected 
via their two-photon decay channel and identified in a standard invariant mass 
analysis from the measured photon momenta. The four photons of an event 
can be arranged in three different combinations to form two 2-photon invariant 
masses (compare Fig. 1). For an acceptable (7, 7r°7r°) event, one of these com- 
binations was required to fullfil the condition, HOMeV < m 77 < 150Me^, 
for both of the 2-photon invariant masses. In addition the mass of the missing 
proton was calculated from the beam energy Eb ea m, target mass m p and the 
energies E n o and momenta p v o of the pions via: 
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Fig. 1. Left: Two photon invariant masses m 7 i 7 2 vs m 7 3 7 4. Right: Missing mass 
Mx — rn p derived from two detected 7r° mesons for incident beam energies Ej < 
400 MeV (symbols with errors: data, histogram: GEANT simulation). 



M\ = ((Ej> + E A ) - (E beam + m p )f - (feo + p^) - (p beam )) 2 (2) 



The resulting distribution is shown in Fig. 1. In case of the reaction ^p — > 
7r°7r°p the missing mass must be equal to the mass of the (undetected) proton 
m„. A Monte Carlo simulation of the 2ir° reaction using GEANT3 (lBru86n 
reproduces the lineshape of the measured data. A cut corresponding to a ±2<r 
width of the simulated lineshape has been applied to select the events of inter- 
est. Random time coincidences between the TAPS detector and the tagging 
spectrometer were su btracted, using events outside the prompt time coinci- 
dence window! Hal96l ). No systematical errors are introduced by this method 



and the statistical errors introduced by the subtraction are included in the 
errors presented in this work. 

The cross section was deduced from the rate of the 27r° events, the number of 
hydrogen atoms per cm 2 , the photon beam flux, the branching ratio of the 7T° 
decay into two photons, and the detector and analysis efficiency. The intensity 
of the photon beam was determined by counting the scattered electrons in the 
tagger focal plane and measuring the loss of photon intensity due to collimation 
with a 100%-efficient BGO detector which was moved into the photon beam at 
reduced intensity. The geometrical detector acceptance and analysis efficiency 
due to cuts and thresholds were obtained using the GEANT3 code and an 
event generator producing distributions of the final state particles according 
to phase space. The acceptance of the detector setup was studied by examining 
independently a grid of the four degrees of freedom for this three body reaction 
(azimuthal symmetry of the reaction was assumed). In a grid of total 1024 bins 
no acceptance holes were found for the beam energy range presented in this 
paper. The average value for the detection efficiency is 1.0%. The systematic 
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Fig. 2. Total cross section for the reaction 7p — » 7r°7r £> (grey squares) at thresh- 
old in comparison with a previous experiment ( WolOOh (open circles) for incident 
energies up to 360 MeV (right) and 425 MeV (left), respectively. T he erro r bars de- 
note the statistical error. The prediction of the ChPT calculation ( Ber96l ) is shown 
(solid c urve) together with its upper limit (dashed curve) and the prediction of Ref. 
i RocOfl (dotted 



curve 



errors are estimated to be 6% and include uncertainties of the beam flux, the 
target length and the efficiency determination. All results presented in this 
work are acceptance corrected and absolutely normalized. 



3 Results and Discussion 



The measured total cross section for the reaction jp — > 7r°7r°p is shown in Fig. 2 
as a function of the incid ent pho ton beam energy. The results are compared 
to a previous experiment ( WolOfj ). The two experiments are consistent within 



the rather large errors of the previous work. 



The prediction of ChPT ( Ber96| ) is plotted up to 40 MeV above the produc- 



tion threshold. The overall shape as well as the absolute magnitude are in 
agreement with the data. Furthermore, the ChPT calculation using the upper 
limit for the Pn(1440) coupling to the s-wave nir channel, can be excluded. In 
the future, the present data might be exploited to provide a better constraint 
on this coupling. Additionally, the cro ss section is compared to the calculation 



with the chiral unitary model (|Roc02j), which especially at threshold predicts 



a smaller cross section. The data show good agreement with both calculations. 

Fig. 3 shows the invariant masses of the 7r°7r° and pir systems in the inci- 
dent energy beam ranges of 330-360 MeV and 360-400 MeV. The pn° mass 
is consistent with a three body phase space distribution, whereas the 7r°7r° 
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Fie. 3. Invariant mass of ir°ir° and 7t°p for different bins of beam energy (full 
squares). The dashed curve shows 3 body phase space. The beam energy range in 
the upper panel is 330-360 MeV and in the lower panel 360-400 MeV. 
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Fig. 4. Angular distributions of the center of mass polar angle O^o (left panel) and 
the angle Vv° (right panel) between the proton and the two-pions in the 7r°7r° rest 
frame (Gottfried Jackson system). The dashed curve shows 3 body phase space. 
The beam energy range in the upper panel is 330-360 MeV and in the lower panel 
360-400 MeV. 



mass deviates slightly already for the energy bin of 330-360 MeV from the 
phase space distributions and shows a tren d towards higher invariant masses. 
The Valencia chiral unitary model ( Roc02| l explains that m^o^o distributions 



skewed to higher invariant masses can arise from the interference of the isospin 
1=0 and 1=2 7r°7r° amplitudes. 

Two angular distributions are depicted in Fig. 4. The polar angle 9 n o of the 
7T° mesons in the overall center of mass frame is consistent with an isotropic 
distribution. The same holds for the angle between the 7r° mesons ifi n o and 
the proton in the frame where the 7r°7r° system is at rest (Gottfried Jackson 
system). Due to the indistinguishability of the two n° mesons, the distribution 
of ipnO shows a symmetry around 90°. The isotropy with respect to the ip n o 
angle in both energy ranges indicates, that the tt° mesons are dominantly in 
an s-wave state. 
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In summary, we have measured the total cross section at threshold for the 
reaction 7p — » 7r°7 r°p. This experiment was motivated as a test of a ChPT 
calculation ( Ber96f ) which shows an unexpectedly high 2tt° rate compared to 



final states including charged pions. This fact is attributed to a dominant con- 
tribution of pion loops which appear in leading (non vanishing) order. This 
prediction is in agreement with our measured data. Furthermore, the upper 
limit quoted for this prediction can be excluded. In future it might be pos- 
sible to exploit the present data to give a better constraint on the Pn(1440) 
to s- wave irn coupling. Secondly, the data are also compared to a predic- 
tion ( Roc02| ). where pion loops are dynamically generated. Especially close to 



threshold, the two predicted cross sections differ significantly. Although the 
present data are of much superior statistical quality than previous measure- 
ments, the precision is still not good enough to discriminate between these 
two models. The observed angular distribution show that the 7r° mesons are 
dominantly emitted in an s-wave state. 
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